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According to the generally accepted scenario, the last giant impact on the Earth formed the Moon and 
initiated the final phase of core formation hy melting the Earth’s mantle. A key goal of geochemistry is to date 
this event, hut different ages have been proposed. Some argue for an early Moon-forming event, approximately 
30 million years (Myr) after the condensation of the first solids in the Solar System \Yin etal.\ |2002| \Jacobsen\ 
\2005^\Taylor et a /.t|2009|l, whereas others claim a date later than 50 Myr (and possibly as late as around 100 My) 
after condensation ( Touboul et al.\ |2007| \Allegre et al.\ |2008| \Halliday\ |2008) . Here we show that a Moon¬ 


forming event at 40 Myr after condensation, or earlier, is ruled out at a 99.9 per cent confidence level. We 
use a large number of N-body simulations to demonstrate a relationship between the time of the last giant 
impact on an Earth-like planet and the amount of mass subsequently added during the era known as Late 
Accretion. As the last giant impact is delayed, the late-accreted mass decreases in a predictable fashion. This 
relationship exists within both the classical sce nario ^O’Brien et al.\\2006\\Raymond et a/.[[2009| l and the Grand 
Tack scenario ^Wakh etal\\2{S\\\\0’Brien cf a/.||2014) of terrestrial planet formation, and it holds across a wide 
range of disc conditions. The concentration of highly siderophile elements (HSEs) in Earth’s mantle constrains 
the mass of chondritic material added to Earth during Late Accretion ^Chyba\\i99\\^ottke et a/.||2010[ l. Using 
HSE abundance measurements pecker et al.\ |2006t \Walker\ [2009) , we determine a Moon-formation age of 
95 ± 32 Myr after the condensation. The possibility exists that some late projectiles were differentiated and left 
an incomplete HSE record in Earth’s mantle. Even in this case, various isotopic constraints strongly suggest 
that the late-accreted mass did not exceed 1 per cent of Earth’s mass, and so the HSE clock still robustly limits 
the timing of the Moon-forming event to significantly later than 40 My after condensation. 

The Moon-forming impact must be the last giant impact experienced by Earth, because both Earth and the Moon 
share an almost identical isotopic composition. This giant impact melted Earth’s mantle and formed the final global 
magma ocean, causing core-mantle differentiation, and it possibly removed a significant portion of Earth’s atmosphere. 
These events can be dated using radiometric chronometers. Unfortunately, the age of the Moon differs substantially 
from one chronometer to the next owing to assumptions in the computation of the so-called model ages. Eor instance, 
Touboul et al. ( 2007|l measured very similar ratios for both Earth and the Moon and, given the differences 


1+90 


in Hf/W ratios estimated at the time for the two bodies, concluded that the Moon-forming event must have been 62 I^q 
M yr after condensation. (Through-out, we use ‘after condensation’ to mean ‘since the birth of the Solar System’; see 
the Supplementary Materials for details.) In this way, the radioactive ^®^Hf would have almost fully decayed into 
beforehand, thus easily accounting for the almost non-existent difference in ratios between the Moon and 
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Earth. However, Konig et al. ( 201 l| l subsequently determined that the HfAV ratios of the Moon and Earth are instead 
identical, voiding this reasoning and leaving the problem of dating the Moon-forming event wide open. 

We approach this problem from a new direction, using a large number of A^-body simulations of the accretion of 
the terrestrial planets from a disk of planetesimals and planetary embryos. The simulations extend across the range of 
well-studied and successful scenarios (for example, typically creating the right number of planets on the correct orbits) 
and fall into two categories. Eirst, classical simulations assume that the disk of terrestrial building blocks extended 
from an inner edge at 0.3-0.7 AU (where 1 AU is the Earth-Sun distance) from the Sun out to 4^.5 AU, just interior 
to Jupiter, and that the giant planets stayed on orbits near their current ones. Our sample of 48 classical simulations 
produces 87 Earth-like planets ^O’Brien et al.\pnQ6\\Raymond et q/.^| 2009| |, which are broadly dehned as hnal bodies 
with masses within a factor of two of Earth’s mass and with orbits between the current orbits of Mercury and Mars. 
In these simulations, the last giant impacts occur between 10 Myr and 150 My after the removal of the solar nebular 
gas, which happened about 3 Myr after condensation. The classical scenario has known shortfalls. The Mars analogs 
are too massive unless the giant planets are assumed to be initially on very eccentric orbits \WetheriU\ 1 1991 [ ^aymond\ 
\et a/.|[2009l l. These large eccentricities, however, cannot be explained in the context of the formation and evolution of 
giant planets in a gas disk. Eurthermore, terrestrial planets accreted in the presence of very eccentric giant planets are 
very low in water content ( Chambers and Cassen 2002[ O’Brien et al. 2006| l. 

The second category consists of simulations in which the terrestrial planets form from a disk with an abrupt outer 
edge at about 1 AU; the inner edge remains the same as in the classical simulations. Simulated solar systems in 
this second category successfully form Mars-like planets ^ansen\ |2009[ \Walsh et al.\ |201 1[ \0’Brien et al.\ |2014| l. 
The truncation of the outer edge can be explained by the early gas-driven inward-then-outward migration of Jupiter 
and Saturn known as the Grand Tack, which then produces a region of greatly depleted surface density between 
the current orbits of Mars and Jupiter \Walsh et al. 201 1[ ). These simulations best reproduce the orbital and mass 
distributions of the terrestrial planets, and they also explain the compositional structure of the asteroid belt ( Walsh 


\et a/.|[20TT] ). Moreover, Earth-like planets accrete volumes of water that are consistent with estimates of the Earth’s 
water content ^O’Brien et al.|[20 T4)l. Previously reported Grand Tack simulations feature last giant impacts occurring 
typically within about 50 Myr after the time of removal of the solar nebular gas Jansen 2009 Walsh et al. 2011 
2014|l. We complement those simulations with new ones (resulting in a total of 211 Grand Tack 
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simulations producing 354 Earth-like planets). 

All Grand Tack simulations produce planetary systems that match the Solar System as well as or better than those 
obtained from classical simulations despite varying many initial disc conditions (see Extended Data Eigure 1), among 
them the initial total mass ratio between the embryo and planetesimal populations (from 1:1 to 8:1, see Supplementary 
Information). If the initial ratio of embryo mass to planetesimal mass is increased, the time of the last giant impact 
also increases (even to about 150 Myr after condensation, see Extended Data Eigure 2) owing to the reduction of the 
well-known effect of dynamical friction—the damping of the eccentricities and inclinations of the larger bodies due to 
gravitational interactions with a swarm of smaller bodies ^O’Brien et a7[[2006) l. Higher eccentricities and inclinations 
of the embryos diminish mutual gravitational focusing, increasing the accretion timescale, and consequently leading 
to later embryo-embryo collisions (that is, giant impacts). 

Dynamical models alone do not indicate whether the Moon-forming impact occurred early (about 30 Myr after 
condensation) or late (about 50-100 Myr after condensation), because the result depends on the initial disc conditions. 
However, we hnd a clear statistical correlation between the time of the Moon-forming impact and the mass subse¬ 
quently accreted, known as the late-accreted mass. This era of Late Accretion includes no giant impacts by dehnition 
and so all of the late-accreted mass comes from the planetesimal population. As shown in Eigure 1, this correlation 
exists across all simulations of both types: classical and Grand Tack. We interpret this correlation by considering that 
the planetesimal population decays over a characteristic time, so that if the last giant impact occurs earlier, then the 
remaining planetesimal population is larger. A larger remaining planetesimal population delivers a larger late accreted 
mass. Strengthening the correlation, a larger initial planetesimal population leads to a shorter timescale for giant im¬ 
pacts owing to enhanced dynamical friction. Eor any given last-giant-impact time, Earth-like planets in the classical 
simulations acquire larger late-accreted masses than those in the Grand Tack simulations (see Eigure 1), because the 
planetesimal population is more dispersed in the classical scenario and therefore decays more slowly. 

The correlation displayed in Eigure 1 can be used as a clock that is independent of radiometric dating systems. The 
late-accreted mass is input into this clock and the time of the last giant impact is read out. A traditional estimate for 
the late-accreted mass can be obtained from the highly siderophile element (HSE) abundances in the Earth’s mantle 
relative to the HSE abundances in chondritic meteorites \Chyba\ 1991[ Bottke et al. 2010| l. HSEs partition strongly 
into iron, and so are transported from the mantle to the core during core formation. In this process, the element 
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ratios are strongly fractionated relative to chondritic proportions [Mann et al. 20121. The HSEs in Earth’s mantle are 
significantly depleted relative to chondritic bodies—a clear consequence of core formation—and yet the remaining 
HSEs are in chondritic or near-chondritic proportions relative to each other pecker et al.||200^ \Walker\ |20091 l. This 
is commonly interpreted as evidence that all or a large portion of the HSEs currently in the mantle were delivered 
by chondritic bodies after the closure of Earth’s core, an accretion phase known as the Late Veneer ( Chou\ 1978] l. 
To account for the observed mantle budget of HSEs, we estimate that a chondritic mass of 4.8 ± 1.6 x 10“^ Mq is 
necessary, where Mq represents an Earth mass. This mass does include contributions from the era known as the Late 
Heavy Bomb ardment. Current mass es timates for this very late (approximately 500 Myr after condensation) accretion 
are 10“"^ M® {Morbidelli et al. 2012a i, which we added to the late-accreted masses of our synthetic Earth-like planets, 
but it only accounts for about 2% of the chondritic mass and therefore does not play a significant role in the analysis 
of the correlation. 

The chondritic mass can only be identical to the late-accreted mass or to the Late Veneer mass if the Moon-forming 
event stripped all of the HSEs from Earth’s mantle or was the last episode of growth for Earth’s core, respectively (as 
is traditionally assumed). However, these conditions are not necessarily true. Consider that some projectiles colliding 
with Earth after the Moon-forming event might have been differentiated, so that their HSEs were contained in their 
cores. If part of these cores had merged with Earth’s core {Albarede et al. 2013| l, then the late-accreted mass would 
clearly be larger than the chondritic mass, because there would be no HSE record of this fraction of the projectile cores 
in Earth’s mantle. Additionally, in this case, given that iron (and therefore HSEs) would have been added to Earth’s 
mantle and its core, the chondritic mass would be larger than the Late Veneer mass, which is geochemically defined 
as the mass accreted to Earth after the core has stopped growing. 

In fact, as explained in detail in the Supplementary Information and Extended Data Eigures 3 and 4, it is unlikely 
that more than 50% of a projectile’s core directly reaches Earth’s core, otherwise geochemical models cannot reproduce 
the tungsten isotope composition of the Earth’s mantle {Rudge et al. 2010| l. Moreover, a large late-accreted mass, 
delivered in only a few objects so as to explain the relative HSE abundances of Earth and the Moon \Bottke etW. 2010| l, 
would have left a detectable isotopic signature on the Earth relative to the Moon, which is not observed ( Wiechert et al.\ 
2001} Zhang et al. 2012| l. Thus, even when considering these more complex possibilities, geochemical evidence 
constrains the late-accreted mass probably not to exceed 0.01 M® (see Supplementary Information). 

For these reasons, we first make the usual assumption that the late accreted mass and the HSE-derived chondritic 
mass are identical. In this case, not a single simulated Earth-like planet with a last giant impact earlier than 48 Myr 
since condensation has a late-accreted mass in agreement with the value estimated from HSEs (see Figure 1). Of 
those forming in less than 48 Myr, only one planet is near the upper 1-cr bounds of the chondritic mass. Only after 
67 Myr since condensation are there Earth-like planets with late-accreted masses consistently within the l-a uncer¬ 
tainty bounds for the chondritic mass. After 126 Myr since condensation, the late-accreted masses of Earth-like planets 
are often significantly below the lower limit set by the HSE measurements. 

We calculate the log-normal mean and standard deviation of the late accreted masses of all Earth-like planets with 
last giant impacts within a range around a chosen time (see Figure 1). We interpret these distributions as a model of 
the likelihood of a specific late accreted mass given a last giant impact time. Given this likelihood model, we compute 
the fraction of Earth-like planets with late accreted masses not exceeding the HSE-derived chondritic mass, also taking 
into account the uncertainties on the latter. Using only the Grand Tack simulations, which provide the best match to the 
terrestrial planet ^Hansen\\2QQ9\\Walsh et al.\\20\\\\O’Bnen et a/.||2014| l, the probability that an Earth-like planet with 
a last giant impact at or before 40 Myr since condensation has an HSE budget consistent with observations is 0.1% or 
less (see Figure 2). Furthermore, if we relax the assumption that the late-accreted mass is equal to the chondritic mass, 
but consider that the former can be up to 0.01 Mq, then a Moon-forming event at or before 40 Myr is still ruled out at 
a 97.5% confidence level or higher. 

Before accepting the proposed HSE clock, a couple caveats need to be discussed. First, our simulations always 
assume perfect accretion—that is, all mass from both colliding bodies ends up in the remaining planet. Accounting for 
imperfect accretion—in which a fraction of the total mass is ejected away from the planet, the true late-accreted mass 
is probably smaller than that inferred from our simulations. Using the characteristics of each impact (for example, 
impact velocity and angle) recorded in our simulations and an algorithm developed from a large suite of numerical 
experiment {Leinhardt and Stewart 2012| l, we estimate the fraction of the projectile mass retained in each impact. 
The effect of imperfect accretion somewhat increases the likelihood of early giant impacts (see Figure 2). However, 
W-body simulations that have incorporated imperfect accretion showed that most of the mass not retained in the giant 
impacts is subsequently re-accreted ( Kokubo and Genda\ 2010 Chambers 20131. Thus, our calculation incorporating 
imperfect accretion should be interpreted as an upper limit to the likelihood for a given last giant impact time. Rein- 
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forcing this analysis, if the HSEs are mostly delivered from the cores of differentiated projectiles as in the scenario of 
Bottke et al. (2010|l, imperfect accretion has a minimal effect on the HSE budget because the material that is lost into 


space comes predominantly from the projectile’s mantle ^sphaug et aZ.|[2006[ ). 

Second, mutually catastrophic collisions between planetesimals may break them into ever smaller pieces until they 
are small enough to be removed by solar radiation before they can be accreted onto planets. If this were a significant 
process, the late-accreted mass would be smaller than estimated from our simulations, which do not include this 
process. However, the size distribution of craters on the lunar highlands suggests that—as in the current asteroid belt— 
most of the planetesimals’ mass was in objects larger than 100 km ^Strom et q/.[| 2005) . The collisional comminution 
for bodies this large is negligible. Moreover, most of the late-accreted mass was probably delivered by Ceres-sized 
(about 1,000 km across) or even larger bodies to explain the relative difference in HSE abundances between Earth and 
the Moon ^Bottke et a(.||201(jl l. 

Considering all of the above, we argue that the clock derived from the correlation between late-accreted mass and 
last-giant-impact age is robust and may be the most reliable way to estimate the age of the Moon-forming event. Given 
the current constraints on the Earth’s late-accreted mass, it establishes that this event occurred significantly later than 
40 Myr and likely at 95 ± 32 Myr after the formation of the first Solar System solids, in agreement with some \Alle gre\ 
et al. 2008t Halliday 20081 (but not all) previous estimates based on radiometric chronometers. Future analyses 


will establish firmer constraints on the late-accreted mass, for instance, by determining the difference in W-isotope 
composition between Earth and the Moon more precisely. From this, we can obtain better limits on the timing of the 
Moon-forming event through our clock. 

Our analysis of numerical simulations can also be used to support or invalidate different accretion scenarios. For 
instance, Albarede et al. ( 2013| l consider the possibility of the impact with Earth of a few projectiles that have a total 
mass of 0.04 M 0 as late as 130 Myr since condensation, but Figure 1 shows that it is extremely unlikely for so much 
mass to be delivered so late. 

A late Moon formation has at least two profound implications. First, it constrains the dynamical conditions of 
the disk from which the planets accreted and the physical properties of the disk material. For instance, a late last 
giant impact implies that most of the mass was in the embryo population rather than the planetesimal population. 
Second, reconciling a late Moon-forming event with radiometric chronometers that suggest the opposite result may 
require challenging fundamental assumptions, such as envisioning a Moon-forming event that did not reset all clocks 
simultaneously and left significant parts of the mantle non-equilibrat ed with the core. This may argue in favor o f 
some of the new scenarios proposed for the Moon-forming collision {Cuk and Stewart 2012 Reufer et al. 2012 1 , 
which distribute impact energy heterogeneously and may leave a significant portion of the Earth’s mantle relatively 
undisturbed. 
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Supplementary Information 


1 Numerical Simulations 

The simulations within the main text come from four published sources ^O’Brien et al.\ |200^ Raymond et al.\ |2009J 
Walsh et al. 2011 O ’Brien et al. \ 20141 and a series of new Grand Tack simulations. They are all N-body simulations, 
which utilize either the Symba ^Duncan et al. 1998|l or the Mercury ( ChambeT^ 1999|l numerical integrator to cal 


culate the gravitational forces between embryos and between embryos and planetesimals (for more details regarding 
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Figure 1: The late-accreted mass relative to each synthetic Earth-like planet’s final mass as a function of the time of 
the last giant impact. Triangles represent Earth-like planets from the first category: classical simulations with Jupiter 
and Saturn near their contemporary orbits ^O’Brien et al.\ |2006[ Raymond et al.\ |2009| l. Circles represent Earth¬ 
like planets from the second category: Grand Tack simulations with a truncated protoplanetary disk ^Walsh et al.\ 
201 1[ O’Brien et al. 20141. The black line resembling a staircase is the moving geometric mean of the late-accreted 
masses in the Grand Tack simulations evaluated at logarithmic time intervals with a spacing parameter of 0.025 and a 
width parameter twice that. The blue region encloses the 1 -ct standard deviation of the late-accreted mass, computed 
assuming that the latter is distributed log-normally about the geometric mean. Always predicting larger late accreted 
masses for each last giant impact time, the dotted staircase is the geometric mean obtained by also considering the 
classical simulations, although the latter do not fit Solar System constraints as well as the Grand Tack simulations do. 
The horizontal dashed line and enclosing darkened region are the best estimate and l-cr uncertainty of the late-accreted 
mass inferred from the HSE abundances in the mantle (chondritic mass): 4.8 ± 1.6 x 10“^ M®. The best estimate 
for the intersection of the correlation and the chondritic mass is 95 ± 32 Myr. The dark and light red regions highlight 
Moon-formation times that are ruled out with 99.9% (40 Myr) and 85% (63 Myr) confidence, respectively. 
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Figure 2; The likelihood that a planet suffering a last giant impact within a specific range of times has a late-accreted 
mass less than or equal to the chondritic mass of 4.8 ± 1.6 x 10“^ M 0 . The probability calculation uses the same bins 
as Figure 1 but only includes the Grand Tack simulations, because they best reproduce the terrestrial planets ^Hansen] 
|2009t \ Walsh et al. \ |201 1[ \0’Brien et al. \ |2014| l. The solid line shows this probability assuming perfect accretion 
and corresponds exactly to the late-accreted masses shown as circles in Figure 1. The lower l-cr limit for the Moon 
formation age is 63 Myr, which corresponds to a 15% probability that an Earth-like planet with a last giant impact 
at that age is characterized by a late-accreted mass equal or smaller than the chondritic mass. The dashed line shows 
the same calculation but for a late-accreted mass less than or equal to 0.01 M 0 , which is an upper limit established 
from a number of elemental and isotopic constraints (see Supplemental Information). The dotted line shows the same 
calculation as for the solid line—that is, using a chondritic mass—^but assuming imperfect accretion during collisions. 
This decreases the late-accreted masses by a variable amount depending on the impact characteristics of the late- 
accreted projectiles onto each planet {Leinhardt and Stewart |2012|l. However this calculation underestimates the 


late-accreted mass because a large fraction of the ejected material would be subsequently re-accreted {Kokubo and 


Genda\ [2010t \Chambers\ |2013|l and projectile core material is less likely to be ejected post-impact ^sphaug et al. 


2006[). Consequently the dotted line overestimates the likelihood that a planet matches the chondritic mass constraint. 


A realistic estimate therefore lies between the solid and dotted curves (probably closer to the former). 
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simulations of terrestrial planet formation see a review by Morbidelli et al. (20T2b|, and for the Grand Tack specifi¬ 
cally, see the Supplementary Information of Walsh et al. ( 201 l| l). The simulations used in this work are broken into 
suites of runs with similar initial conditions. A suite of simulations consists typically of 10, but anywhere between 
8 and 20, independent runs with similar initial protoplanetary discs. In the main text, we identified two categories of 
simulations that all of the suites fall into: (1) classical and (2) Grand Tack. 

The classical model of terrestrial planet formation assumes that Jupiter and Saturn formed on or near their current 
orbits. They did not undergo any large-scale migration during the few million year gaseous disc phase and were fully- 
formed during the last phases of terrestrial accretion. The N-body simulations begin after the gas has been dispersed 
(i.e. ~ 3 My after the first solids in the Solar System ( Haisch et aI.||200T] l). The simulations used in the main text are 
taken from O’Brien et al. ( |2006| l and Raymond et al. (2009 i. The initial conditions of these simulations include solid 
bodies from an inner edge at 0.3-0. 7 AU out to 4^.5 AU. They include a few tens to ^ 100 planetary embryos (Moon- 
to Mars-mass bodies) and ^ 1000-2000 planetesimals (Ceres-mass bodies) with an equal mass in each population. 
The orbits of the giant planets are varied in different simulations, but Jupiter is always close to its current location. 
The most important parameter that was changed in these simulations was the giant planets’ eccentricities, which were 
tested between 0-0.1 . Extended Data Figure 1 shows as open triangles the outcome of these simulations including the 
mass-orbit distribution, the timing of the last giant impact on Earth, the late accreted mass, and two orbital structure 
statistics. 

We note that simulations using the classical model have consistently failed to reproduce the terrestrial planets’ 
mass-orbit distribution \Wetheriii\ |1991[ \Chamber^ |2001t \0’Brien et al.\ |2006[ {Raymond et al.\ |2006[ |2009[ \Mor-\ 
\ishima et al.\ |20T0|. Specifically, simulations that are consistent with the current larger-scale picture of the Solar 
System’s evolution (i.e., including a giant planet i nstability (|LevAo« et al. 201 l| l) produce Mars analogs that are an 
order of magnitude too large {O’Brien et al. 2006[ Raymond et al. 2009| l (shown as open triangles in panel a of Ex¬ 
tended Data Figure 1). Simulations with giant planets on more eccentric orbits can produce improved Mars analogs, 
however even these improved Mars analogs are often still a factor of two too large at 1.5 AU (also shown as open tri¬ 
angles in panel a of Extended Data Figure 1). Furthermore, classical simulations that produce improved Mars analogs 


are inconsistent with the late planetesimal-driven migration of the giant planets ( Raymond et al. 2009 Levison et al. 


2011 


et al. 


and also tend to form dry terrestrial planets on overly-excited orbits ^Chambers and Cassen 2002 Raymond 


[20091 1. 


These problems served as the main motivation for the Grand Tack scenario, the second category of simulations. 
Regardless of the initial conditions, which vary as described in the next paragraph, the Grand Tack simulations always 
proceed as follows. Similar to the classical model, each simulated terrestrial protoplanetary disc starts with a disc of 
~ 100 embryos and ^ 1000-2000 planetesimals. The Grand Tack simulations begin slightly earlier than the classical 
simulations at 0.6 My before gas dispersal (i.e. ~ 2.4 My after the first solids in the Solar System, if the disk disappears 
after 3 My {Haisch et al. 200 Ij l). \Walsh et~ar. ( 201 l| l did an extensive exploration of possible giant planet migration 
evolutions (timescales, initial locations, role of Uranus and Neptune, etc.) and concluded that the sculpting of the 
inner protoplanetary disc is remarkably insensitive to these details. The key factor is the location of Jupiter’s “tack,” 
which occurs when Jupiter and Saturn enter into resonance and change the direction of their migration from inward to 
outward. Therefore, we use the simplest scheme tested h y \Walsh et al. (2011 1 namely, “Saturn’s core growing in the 
2:3 resonance with Jupiter” (see Walsh et al Walsh et al. ( 2011| l supplementary material). In this scheme, during the 
first 0.1 My, Jupiter and Saturn migrate inward from 3.5 and 4.5 AU to 1.5 and 2 AU, respectively, while Saturn’s mass 
grows linearly from ~ 10 M® to its current mass. This gas-driven migration truncates the embryo and planetesimal 
disc at ~ 1 AU near the 3:2 mean motion resonance with Jupiter. When Saturn reaches a mass close to its final mass, 
the migration physics changes and Jupiter and Saturn migrate outwards to 5.25 and 7 AU, respectively, in 0.5 My. 
At this point the gas is removed and all giant planet migration stops. These final locations for Jupiter and Saturn 
are appropriate initial conditions for the late giant planet instability ^orbidelli et al. |2007| l. Then we evolve the 
simulation for another 150 My. At the end of this period, there are ~ 4 terrestrial planets on stable orbits. 

(|201 1|| and 16 additional 


For the second category of the simulations, we use the 8 simulations from Walsh et al. 


simulations fromjO ’Brien et q/^(|2014|i, which are statistical variants of the\Walsh et al. (20111 runs with the addition of 


outer planetesimals scattered inward by the outward migration of the giant planets. We also use 187 new simulations. 
The Grand Tack simulations can be grouped into 27 suites. The initial protoplanetary disc conditions are very similar 
to \Walsh et a/.| ( |20TT| l with the exception of two parameters: the ratio of the initial amount of mass in embryos to the 
amount of mass in planetesimals and the initial mass of each individual embryo. We also explored the consequences 
of an initial distribution of embryos with mass increasing with semi-major axis. 

Regardless of the varied parameters, the mass-orbit distributions are much better than that obtained in the classical 
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HSE 

Mantle 

concentration. Cm 

Chondritic 
concentration, Cc 

Mantle to 
chondritic ratio, R 

Re 

3.5 ±0.6 X 10-1° 

5.5 ± 1.5 X 10"® 

6.3 ±2.0 X 10-° 

Os 

3.9 ±0.5 X 10"° 

6.5 ± 1.8 X 10-1’ 

6.0 ± 1.8 X 10-° 

Ir 

3.5 ±0.4 X 10"° 

5.9 ± 1.9 X 10-1’ 

5.9 ±2.0 X 10-° 

Ru 

7.0 ±0.9 X 10"° 

8.7 ±2.5 X 10-1’ 

8.0 ±2.5 X 10-° 

Pt 

7.6 ± 1.3 X 10"° 

1.2 ±0.3 X 10-° 

6.5 ±2.1 X 10-° 

Pd 

7.1 ± 1.3 X 10"° 

7.3 ±2.2 X 10-1’ 

9.7 ±3.4 X 10-° 


simulations (see panel a of Extended Data Figure 1). In panel b, we show the primary result of the main text that the 
relative late accreted mass is correlated with the time of the last giant impact, and that the Grand Tack scenario gives 
on average smaller late accreted masses at any given last giant impact time. We note that Grand Tack simulations are 
consistently more concentrated than the classical simulations independent of the relative late accreted mass (panel c), 
and are generally as or less dynamically excited than classical simulations (panel d). Nevertheless, those simulations 
that are consistent with estimates of the late accreted mass (gray horizontal band in panels b-d) produce terrestrial 
planet systems, which are in general too excited and less concentrated than the real terrestrial planets. This is likely 
due to the assumption of perfect accretion during collisions, which means that the role of impact ejecta is not accounted 
for in the simulations. In fact, it has been shown that impact ejecta can significantly reduce the final excitation of the 
terrestrial planets ( | Chamber^ |20 13 [ ). This would reduce the angular momentum deficit and increase the concentration 
parameter, since dynamically colder systems are capable of maintaining a more compact configuration as evidenced 
by the loose correlation between Sd and Sc- Therefore, one should not interpret panels c and d as definitive evidence 
that the relative late accreted mass should be high, since it’s likely that implementing imperfect accretion will change 
both Sd and Sc metrics for any given late accreted mass. 

Of all the parameters that we changed in the Grand Tack simulations, we found that only the ratio of the total mass 
of embryos to planetesimals affects the last giant impact time and the late accreted mass (see Extended Data Figure 2). 
It is clear that by varying the initial disc conditions it is possible to obtain a wide range of last giant impact times and 
even obtain a wide range of last giant impact times within the same simulation suite. But, regardless of the simulation 
suite, each last giant impact date corresponds to only a tight range of late accreted masses. 

Naively, it is surprising that the aforementioned range is so tight, since one might think that by changing the initial 
mass in planetesimals in a simulation, one could get a similar last giant impact time with a different late accreted 
mass. However, this is not the case because the mass in the planetesimal population governs dynamical friction, which 
determines when the giant impacts occur. Thus increasing the planetesimal population makes the giant impacts occur 
earlier and decreasing it makes them occur later. In other words, changing the planetesimal mass makes the final 
results evolve along a diagonal line (top-left to bottom-right) in Figure 1. This is the reason why for a given last giant 
impact date, the spread in late accreted mass values is relatively narrow. 

The difference between the correlations given by the classical and Grand Tack simulations is due to a fundamental 
change in the structure of the disk. During the Grand Tack, Jupiter rapidly clears the outer terrestrial region (i.e. 
asteroid belt and Mars region). Instead, in the classical simulations, this region is de-populated slowly due to resonant 
interactions with Jupiter. This leads to a slower decay of the planetesimal population and hence a larger mass delivered 
for a given last giant impact time. We notice however that the dispersions of late accreted mass for a given last giant 
impact time are comparable for the two scenarios. 

2 Considerations on Late Accretion 

In the main text of this paper we showed the existence of a correlation between the timing of the Moon-forming event 
and the amount of material accreted after this event, known as the late accreted mass. In order to date the Moon¬ 
forming event using this correlation, we need a reliable constraint on the late accreted mass. Here we discuss various 
constraints, which support the usual assumption that the late accreted mass is comparable to the usual estimate of the 
late veneer mass, i.e. ~ 5 x 10“^ M 0 . Nevertheless, the possibility exists that late accreted mass could have been 
somewhat larger (but unlikely exceeding 1% of an Earth mass as we will see below). We proceed in steps below, first 
presenting a definition of Fate Accretion, then introducing a set of compositional and isotopic constraints. 
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2.1 Introduction to Late Accretion: 

After the Moon-forming impact, the Earth continued to grow. In the context of our simulations, we define the late 
accreted mass Mla as simply the planetesimals accreted by the target Earth-like planet after the last giant impact. A 
giant impact is a collision between the Earth and a surviving embryo (i.e. oligarch) whereas planetesimals represent 
the surviving small body population. Embryos are several thousand kilometers in size (i.e. Moon to Mars sized) and 
planetesimals range up to, possibly, a couple thousand kilometers in size (Ceres, the largest planetesimal remaining in 
the asteroid belt has a diameter of 900 km). The bimodal mass distribution of embryos and planetesimals is a direct 
outcome of the oligarchic growth regime of terrestrial planet formation ^Kokubo and Ida 1998|l. 
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20121. Once the Earth’s mantle solidified in the aftermath 




of the Moon-forming event, the planetesimal bombardment and the remnant radioactive elements were probably not 
enough of an energy source to re-melt the planet completely. Thus the fate of iron and siderophile elements delivered 
by late projectiles is determined by the nature of the projectile and the characteristics of its impact. If the largest 
planetesimals are differentiated, their cores may partially merge with the Earth’s core ^Olson and Weeraratne\^QQ^ . 
Thus, the termination of terrestrial core formation is ambiguous, i.e. probably not an abrupt or global event. This 
confuses the link between between Late Accretion and Late Veneer, the latter being a geochemically defined entity, 
which refers to the mass accreted by the Earth after core closure ^Dauphas and Marty\\2QQ2\\Maier et al.||2009HMann 

[iFHri[20T2l l. 

In order to better constrain the late accreted mass, we assume that all planetesimals, differentiated or not, are bulk 
chondritic. An undifferentiated planetesimal delivers a mass of chondritic material equal to the planetesimal’s mass. 
Eor a differentiated planetesimal, a fraction X of the mass of the core may emulsify. In order to explain the excess 
of HSEs in the Earth relative to the Moon, the emulsified droplets have to be permanently integrated into the Earth’s 
mantle ( Bottke et al. 2010| l, presumably because the oxidizing conditions in the mantle prevent the iron from reaching 
the core. The rest 1 — X of the mass of the projectile’s core sinks as a competent blob through the mantle and merges 
with the Earth’s core. The entire mantle of the projectile goes into the Earth’s mantle. Eor convenience, we can split 
the projectile’s mantle into two parts as well with proportions X and 1 — X. In this way, it is clear that a fraction X of 
the total mass of the planetesimal is delivered as bulk chondritic material to the Earth’s mantle. The remaining fraction 
1 — X of the projectile’s mantle provides an achondritic contribution to the Earth’s mantle. 

This simplification requires us to assume that during Late Accretion no significant mass is delivered by achondritic 
bodies, i.e. bodies without the metal counterpart that would restore a bulk chondritic composition. These are probably 
ejecta from the mantles of differentiated embryos or the Earth itself, launched into space during the Moon-forming 
event ^eufer et al 2012| l or previous hit-and-run and erosive giant impacts ^sphaug et ah 200()\^hambers 2013| l. 
This achondritic mass is difficult to constrain using any elemental or isotopic system. Eortunately, we do not need 
to worry about these fragments. Our simulations address only the chondritic component of late accreted material 
because the ejection (and potential re-accretion) of mantle material in hit-and-run collisions is not taken into account. 
This means that we can apply our correlation between the timing of the Moon-forming event and the amount of late 
accreted material considering solely the mass delivered by bodies of bulk chondritic properties. 

2.2 Highly siderophile element constraints: 


Planetary differentiation within a melted mantle is very efficient {Ruble et al.\ 2004 1 , therefore each giant impact 


including the Moon-forming event, sequesters into the core a large fraction of the highly siderophile elements (HSEs) 
delivered prior to the impact. Any remaining HSEs in a melted portion of the mantle are strongly fractionated relative 
to chondritic by the differentiation process as well as being significantly depleted ^Mann ef ar||2012[ ). Once the mantle 
has solidified after a giant impact, it is much more difficult for highly siderophile elements to enter the core unless the 
mantle is locally melted during the impact of a large planetesimal. In this case, it is possible for the projectile’s core 
to entrain HSEs already resident in the Earth’s mantle into the Earth’s core. However HSEs in the Earth’s mantle are 
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Mantle concentration 

Chondritic concentration 

Mantle to chondritic 

HSE 

Cm [ng/g] 

Ce [ng/g] 

ratio R 

Re 

0.35 ±0.06 

55.3 ± 15.0 

6.3 ±2.0 X 10-3 

Os 

3.9 ±0.5 

653 ± 180 

6.0 ± 1.8 X 10-3 

Ir 

3.5 ±0.4 

592 ± 185 

5.9 ±2.0 X 10-3 

Ru 

7.0 ±0.9 

872 ± 253 

8.0 ±2.5 X 10-3 

Pt 

7.6 ± 1.3 

1174 ±334 

6.5 ±2.1 X 10-3 

Pd 

7.1 ± 1.3 

732 ± 222 

9.7 ±3.4 X 10-3 


Table 1: Chondritic concentrations for different HSEs (column 1) in the Earth’s mantle ^Becker et a/.||2006l l (column 
2) and chondritic meteorites ^Walker\\2QQ9\ (column 3) and the ratio of the two (column 4). 


observed to have chondritic proportions to one another [Becker et al. 2006 Fischer-Godde et ar\ 201 l| l (see Table[2. 
This implies that the vast majority of the HSEs today in the terrestrial mantle did not fractionate during metal-silicate 
segregation. Reaffirming this hypothesis, the ^®^Os/^®®Os isotope ratio (one of the HSEs) is chondritic for both the 
Earth and the Moon ^eisel et a/.[|1996l|iyafer et a/.|[2004{|lTQ/^^|20091 l. 

The concentrations and associated l-cr uncertainties of 6 highly siderophile elements [Becker et al. 2006] l in the 
Earth’s mantle are reported in Table[2 We also include the ‘chondritic mean’ and its uncertainty for each of the HSEs, 
that we determined from the mean, minimum and maximum values for ordinary, carbonaceous and enstatite chondrites 
reported in \Walker\ ( |2009[ ). We define the mean chondritic concentrations as the average of the mean values for the 
three chondrite groups. We estimate the uncertainty by assuming that half the difference between the maximum and 
the minimum value for each chondrite group is 3 -(t, then average the estimated l-cr values of each chondrite group to 
obtain the uncertainty of the chondritic concentrations. 

The chondritic mass needed to account for the HSE terrestrial mantle budget is Mq = R x Mm where R is the 
ratio of mantle to chondritic abundances of each HSE and Mm = 0.675 Mq is the mass of the Earth’s mantle. Each 
HSE provides its own chondritic mass estimate, and the l-cr uncertainty is determined by standard error propagation. 
Averaging them together, we obtain a chondritic mass of Me = 4.8 ± 1.6 x 10“^ Mq. 

This chondritic mass corresponds directly to the amount of chondritic material necessary to account for the mea¬ 
sured HSE abundances in the Earth’s mantle relative to chondritic. We purposely avoid identifying this chondritic 
mass with either the late veneer mass or the late accreted mass. By definition, the late veneer mass was delivered after 
core formation ended, and the late accreted mass was delivered after the Moon-forming event. Only under a specific 
set of assumptions are either of those labels accurate for the chondritic mass, and so we will be careful delineating 
when different sets of assumptions are being made. 

The simplest set of assumptions we can make is to assume that either all late projectiles were undifferentiated or 
their cores fully dissolved in the Earth’s mantle X = 1 with no metal percolating to the core (pre sumably because of 
the oxidized state of the mantle via reaction with [Frost et al. 2004 Frost andMcCammon 2008 \ Ee^+ or through the 
oxidizing effect of water ^Sharp ef a/.||2013l l), or some combination of the two. In addition, we assume that the Moon¬ 
forming event cleared the entire terrestrial mantle of its previously existing HSEs. Only if all of these assumptions are 
true, does the late accreted mass Mla coincide with the measured chondritic mass based on HSEs: 


Mla = Mc= 4.8 ± 1.6 X 10"^ M, 


In this scenario, no late projectile cores merge with the Earth’s core and no mantle HSEs pre-date the Moon-forming 
impact. This is the assumption most often made in the literature \Bottke et al.\ |2010t |Ma«« et al.\ |2012[ \Wang ari^ 
\Becker\ |2013| l and is the nominal case we present in the main text. Below, we expand this model to a more general 
case that allows the late accreted mass, the late veneer mass and the chondritic mass to be distinct. 

Before we discuss the fate of the projectile cores, we need to assess what portion of Late Accretion was delivered 
in large bodies with cores. We do this by considering that during Late Accretion, HSEs accumulate not only on the 
Earth but also on the Moon ^Chyba\\\99\\ . However, the lunar chondritic mass ( |Day et a7r||2007[|£)Qy and Walker\ 
20111 is me ~ 1.8 x 10“® Mq; it is computed in the same way, from a bulk HSE abundance ^ 0.00015 x chondritic). 


It should be kept in mind that the concentration value for the Moon is more uncertain than for the Earth because it 
has been obtained from indirect sampling of the lunar mantle (through igneous rocks such as basalts). The ratio of 
lunar to Earth chondritic mass is Me/me « 2700, whereas the ratio of accreted material according to the ratio of 
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gravitational cross sections is cr « 13.5 +1610/ (5.6 + V^), where 14o is the relative velocity of the projectile before 
any acceleration from the gravitational potential of the target. The expected flux ratio varies from a « 300 for bodies 
with essentially no relative velocity with the Earth-Moon system (i.e. dust), to cr « 52 for asteroids with 140 ~ 6 km 
s“^ and (T « 15 for comets with V/o ~ 40 km s“^. There is a discrepant factor of 9, 50 or 200 between the measured 
and the expected ratios of chondritic masses delivered to the Earth and the Moon, depending on the projectile relative 
velocity 14o = 0, 6, or 40 km s“^, respectively. 

To reconcile this discrepancy with the assumption that the chondritic mass is entirely delivered during Late Ac¬ 
cretion, \Bottke et al.\ ( |2010| l postulated that most of the mass accreted by the Earth was delivered by a few large 
planetesimals. Since the ratio of gravitational cross-sections only dictate a statistical long-term average of the impact 
rates, stochasticity and small number statistics can explain why the few largest planetesimals delivered mass only to 
the Earth and not to the Moon {Raymondet al. 2013| l. Eor instance, only one 2500 km diameter planetesimal is needed 
to deliver 85% of the chondritic mass assuming a density of 3 g cm“^. If the largest planetesimals were only 1500 km 
in diameter, then five such bodies could deliver 93% of the chondritic mass. We quantify this hypothesis by dividing 
projectiles into two size bins: large Ml and small Ms- Thus, Mla = Ml + Ms- Large projectiles are rare and 
thus they are likely to be accreted only by the Earth and not by the Moon, while small projectiles are numerous and 
therefore they are delivered to the Earth and the Moon proportionately according to their respective cross-sections a. 


Willbold et al. (20111 also finds evidence, using measurements of the Earth’s tungsten anomaly, for a ‘patchy’ Late 


Accretion consistent with a few large impacts rather than many smaller impacts. 

The \Bottke et aL\ ( |2010| l scenario is simple and matches the data, but its hypothesis rests principally on two as¬ 
sumptions, both of which are responsible for the conclusion that the late accreted mass is equivalent to the chondritic 
mass derived from the HSE abundances of the Earth’s mantle relative to chondritic meteorites. Eor a more accurate 
assessment of the relationship between the late accreted mass and the chondritic mass, we identify each assumption 
and adjust the model accordingly. 

Eirst, we can challenge the assumption that all of the HSEs were delivered late (i.e. after the Moon-forming 
impact). The survival of mantle HSEs during the Moon-forming event is a possibility to be taken seriously. On 
the modeling side, new high-resolution smooth-particle hydrodynamic (SPH) simulations of the Moon-forming even t 


show that the opposite hemisphere of the Earth can remain unaffected by the giant collision ( Cuk and Stewart 


20121 . 


While HSEs in the melted hemisphere are transported into the core and are strongly fractionated relative to chondritic, 
they are also strongly depleted relative to the other hemisphere, where the HSEs retain their chondritic relative propor¬ 
tions. Thus, in this simple scenario, in the end the bulk mantle of the Earth contains HSEs in approximate chondritic 
relative proportions and depleted roughly by a factor of 2 relative to their pre-Moon-forming-impact abundances. On 
the geochemical side, Touboul et al. (2012 1 discuss Kostomuksha komatiites, for which combined ^®^W,^®®487Qg^ 


i42,i43]vjj isotopic data indicate that their mantle source underwent metal-silicate fractionation well before 30 My of 
Solar System history and therefore remained unaffected by the Moon-forming event, assuming that the latter occurred 
not earlier than this date. This supports the scenario described above. Thus, the HSE budget may be a combination of 
surviving and late delivered HSEs, and so we quantify the relationship between the chondritic mass Mq corresponding 
to the mantle budget of HSEs and the chondritic component of the late accreted mass Mqla by introducing a param¬ 
eter Y, which represents the fraction of HSE’s delivered during Late Accretion. This parameter can range between 
Y = 1, if nothing pre-dates the Moon-forming event, to L = 0, if everything does. Thus the chondritic component of 
the late accreted mass is related to the chondritic mass determined from the HSEs: Mqla = YMq. 

Second, we can challenge the assumption that that the projectile cores merged entirely with the Earth’s mantle. 
This requires that the core is emulsified or at least broken into incoherent pieces that interact efficiently with the 
mantle l^ubie et al. 2003 i and oxidized via reaction with Ee®+ {Frost et al. 2004|| Frost and McCammon 2008 


Bottke et al. 


2010|l or alternatively with water ( Sharp et al. |2013|l to prevent the transportation of iron and siderophile 


elements into the Earth’s core. Recent results suggest that it is possible for part of the projectile’s core to remain intact 
and merge directly with the Earth’s core. Using smooth-particle hydrodynamics simulations, Robin Canup (personal 
communication) finds that only up to half of the largest projectile cores can merge with the Earth’s core. On the 
other hand, physical experiments ( Olson and Weeraratne\ 2008[ Deguen et al. 201 1)| and numerical simulations {Dahl 


and Stevenson 2010[ Samuel 20121 of metal sinking through the mantle conclude that projectile cores with radii 


much smaller than the depth of the magma ocean emulsify efficiently. Erom these studies, it is unclear whether larger 
planetesimal cores may survive the descent towards the core-mantle boundary. As we already said, we parameterize 
this uncertainty by asserting that an average fraction X of each core dissolves in the mantle of the Earth, while the 
remaining fraction 1 — X merges into the Earth’s core. The X parameter can range between 1 (complete assimilation 
of the projectile’s core into the Earth’s mantle) and 0 (the entire projectile’s core sinking onto the Earth core). We 
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stress that the fraction X represents the mass weighted average for all the late accreted large planetesimals, not each 
individual one: 


N 


X = 


Ml 




where N is the total number of large planetesimals, Mi and Xi are their individual masses and values of X. In this 
average, if large undifferentiated planetesimals, for instance very oxidized bodies, contribute to Late Accretion, then 
they count as if they have Xi = 1. 

We had already divided the late accreted mass into components that were delivered either by small and large bodies: 
Mla — Ml + Ms- Reminding ourselves that we have assumed that all bodies are chondritic in bulk, small bodies 
thus contribute their entire mass to the chondritic component of the late accreted mass. Large and hence differentiated 
bodies only deliver on average a fraction X of their mass to the chondritic component of the late accreted mass. Thus, 
the chondritic component of the late accreted mass is: 


Mcla = XMl + Ms = Y Mo¬ 


using this set of relationships and the expression of the late accreted mass as the sum of two size components: Mla = 
Ml + Ms, we can re-express the late accreted mass as a function of the unknown parameters X and Y: 


1 ( TTXc^ 

Mla = —{YMc - crmc) + crmc = ( 1 - 


1-X 

Y 


where Me = 4.8 x 10“^ is the chondritic mass in the Earth’s mantle, me = 1.8 x is the lunar late accreted 

mass, a is the cross-section ratio between the Earth and the Moon, and Ms = erme is the mass delivered by small 
projectiles to the Earth. Because a is much smaller than the ratio Me/me the following analysis is not sensitive to 
the choice of a. 

Extended Data Eigure 3 shows the late accreted mass Mla function of X and Y given the teiTestrial Me and 
lunar chondritic masses me, and a gravitational cross-section ratio tr = 52 appropriate for planetesimals with semi¬ 
major axes in the terrestrial region. As the parameters X and Y change, the mass delivered in large planetesimals 
must adjust to continue to match the HSE constraints. Immediately, we see that a late accreted mass of 5 x 10-3 Me 
approximately divides the parameter space in half, mnning diagonally from near the Bottke et al. ( |2010 l scenario, 
which is exactly at Ai = 1 and Y = 1, to a scenario where all of the HSEs pre-date the Moon-forming impact and the 
core of each projectile accretes entirely into the Earth’s core X = Y = 0. Scenarios below this contour in the hgure 
require smaller late accreted masses than the \Bottke et a/.|p010| l scenario, while those above require more. Most of 
the parameter space has a late accreted mass Mla < 0.01 Mq. 

2.3 Constraining X: 

Up to now, X has been treated as a free parameter; as X approaches zero, the late accreted mass can be arbitrarily 
large. We now try to constrain X by addressing what fraction of the projectile’s core emulsihes after impact with 
the Earth. The emulsihcation of the projectile’s core is a mechanical (hydrodynamical) process that does not depend 
on the chemical properties of the Earth, which instead determine the fate of the metallic droplets. Eor this reason, 
we can use works that constrain the fraction of the projectile emulsihed in an impact \Rudge et al.\ 2010[ Dahl 
and Stevenson |2010[ Deguen et al. 


2011 Samuel 2012} regardless of what these work envision happening to the 


droplets. The droplets can percolate into the Earth’s core if the mantle is sufficiently reduced as in the early phases of 
accretion ^ubie et a/.||201l) or be oxidized and remain in the mantle as necessitated to explain the HSE abundances 
and their chondritic relative proportions in the Earth’s mantle [Bottke et al. 2010| l. Physical experiments [Olson and 
\Weeraratn^ |2008[ \Deguen et al.\ |2011| l and numerical simulations ( Dahl and Stevenso^ |2010[ \Samuel\ |2012| l show 
that emulsihcation is important until the projectile core approaches the thickness of the mantle; however, it is difficult 
to derive a lower bound on X given the problems in modeling the emulsihcation process which involves length scales 
ranging from centimeters to hundreds of kilometers or more. Thus, we think that the best constraints on X can be 
derived from geochemical analysis ^udge et al 20101. 

Specihcally, we apply the model of Rudge et al. (20^ for the evolution of the mantle’s tungsten isotopic anomaly 
during the accretion history. The coefficient k in Rudge et al. ( |2010 1 determines the fraction of the projectile cores 
that emulsify in the Earth’s mantle like our coefficient X, so they have the same meaning for our purposes. 

\Rudge et a/.| ( |2010'] l showed that isotopic constraints from the Hf-W and U-Pb systems bound the accretion history 
of the Earth. If growth is too fast or too slow, too much or too little radiogenic tungsten and lead accumulate in the 
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Earth’s mantle. The bounds presented in Rudge et^ ( 2010| l assume constant partitioning coefficients, which in reality 
are a function of temperature, pressure and oxygen fugacity. The oxygen fugacity can change significantly during the 
early accretion history ^Wade and Wood 2005 1 Rubie et ar\ |201 1[ ) while the temperature and pressure are required to 
lie on the peridotite liquidus. However, Rudge et al. ( 2010| l discovered that isotopic evolution models with and without 
evolving partition coefficients do not differ substantially. The only difference they found when using evolving partition 
coefficients “is that accretion needs to be slightly more protracted to match the same Hf/W and U/Pb observations. 
The requirement of more protracted accretion arises because both W and Pb are more siderophile during the early 
accretion than the Late Accretion, which causes a bias towards younger ages.” 

Panel a in Extended Data Eigure 4 shows the X > 0.53 bound determined following Rudge et al. ( 2010[) from 
Hf-W constraints (U-Pb does not constrain the equilibration factor). This bound is not identical to that stated in^udge 


et al. P010| l (0.36), since we have used the updated terrestrial HfAV measurements of |A'd«/g et al. ( 201 l| l. Examining 
panel a, we see that all combinations of AT > 0.53 and Y produce late accreted masses M^a < 8.9 x 10“^ Mq. 


2.4 Isotopic Constraints 


Additional constraints on the late accreted mass of the Earth come from the fact that the Earth and the Moon have 
very similar isotopic contributions. If a large mass of a specific composition had been added to the Earth after the 
Moon forming event, then the Earth and the Moon would be more different than observed. Here we first examine the 
isotopic systems of oxygen and titanium, which show that the terrestrial late accreted mass is unlikely to exceed 1% of 
an Earth mass regardless of the dominant projectile composition with the possible exception of enstatite chondrites. A 
differentiated enstatite chondrite however should have a strong radiogenic tungsten signature given its reduced nature. 
Therefore, we next examine constraints provided by the lunar-terrestrial tungsten anomaly, which is the difference 
between the tungsten isotopic composition of the Moon and Earth. 

Oxygen and titanium isotopic systems: If most of the late accreted mass was delivered by one or a few large differenti¬ 
ated planetesimals, as required to explain the difference in HSE abundances of the Earth and Moon, its chemical nature 
could be easily dominated by a single parent body composition. Different isotopic systems can constrain contributions 
to the late accreted mass for projectile compositions analog to most meteorite types. 

\Wiechert et a/.| ( |200T] l and \Spicuzza et al.\\2001\ found that the Earth and the Moon lay on the same oxygen isotope 
fractionation line A^^Omooh -- A’^’^OEm-th = 0.008 ± 0.01%o where A^’^Osampie = <5^^0sampie — 0.5245 x (5^®Osampie 
and (5*Osampie = 10 ^ X [(*0/^®0)sampie/(*0/^®0)standaid — l] ■ All reported uncertainties are 1-cr. Recently, Herwartz 
and colleagues announced at the Royal Society/Kavli Institute meeting on the Origin of the Moon - challenges and 
prospects (25-26 September, 2013) at Chicheley Hall, Buckinghamshire, U.K. that they resolve a difference between 
the two bodies; A^^Omooh ~ A^^OEauh = 0.012%o. If this measurement is correct and if we assume that the Earth 
and Moon equilibrated oxygen isotopes at the time of the Moon-forming event, knowing that the Earth received more 
late accreted mass than the Moon, we conclude that the dominant nature of the projectiles during Late Accretion must 
have had a carbonaceous chondrite (excluding Cl), enstatite or HED composition, because these are the only meteoritic 
compositions that are situated below or on the terrestrial fractionation line in the oxygen-isotope diagram. 

Then we can constrain how much late accreted mass could be added to the Earth from carbonaceous CV, CO, CK, 
CM, CR and CH meteoritic compositions using their measured oxygen isotope compositions ^Clayton and Mayeda 


19991: Ai^Ocv = -3.3 ± 0.6%o, Ai^Oco = -4.3 ± 0.3%o, A^^Ock = -4.2 ± 0.4%o, AI^Ocm = -2.3 ± l.l%o. 


A^^Ocr = —1.6 ± 0.5%o and A^^Och = -1.5 ± 0.2%o. A late accreted mass dominated by each composition is 
limited to be less than 0.4 ± 0.3%, 0.3 ± 0.2%, 0.3 ± 0.2%, 0.5 ± 0.5%, 0.8 ± 0.7% and 0.8 ± 0.7% of an M®, 
respectively. Eor HEDs and enstatites, oxygen isotopes are not very useful, because these meteorite groups are much 
closer to the terrestrial fractionation line. 

Euither constra ints are provided b y the relative titanium isotope composition of the Earth and Moon. The Moon- 
Earth difference is (Zfiangef a/. 20121 e®°TiMoon—e^°TiEarth = —0.04±0.02 where e®°Tisampie = lO'^x [(^°Ti/^'^Ti)sampie/(®°Ti/"‘’^Ti)stai 
1]. HEDs and, to a lesser extent, enstatites, can be excl uded as the domina nt component of Late Accretion be- 

1201211 (e^OEWstatite = -0.23 ± 0.09 and 


cause they have negative titanium isotope compositions [Zhang et al. 


e^^TinED = —1.24 ± 0.03) and any large addition of this composition would result in a positive e^'^Ti measure¬ 


ment difference between the Moon and the Earth. A source of uncertainty in this analysis appears from the possibility 
that the Moon-Earth difference is positive, since the 2-a uncertainty encompasses the zero value. Considering this 
possibility and assuming that e®°TiMoon — e^°TiEaiih = 0.01 (a 2.5 -(t deviation), we can still constrain how much late 
accreted mass could be added to the Earth from HED meteoritic compositions to 0.8% of an M®, however we cannot 
place such a strict constraint on an enstatite chondrite composition. 

In summary, oxygen and titanium isotopic constraints limit the dominant component of the late accreted mass to 
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less than 1% of an Earth mass for most meteoritic compositions. The only exception may be an enstatite chondritic 
composition, if the Earth-Moon difference is e^^TiMoon — e^^TiEanh > 0.002, which is barely exterior to the upper 
2-a uncertainty of the current measurements, then a considerable mass with enstatite composition could be added to 
the Earth. However, we can probably still exclude differentiated enstatite chondritic bodies from dominating Late 
Accretion. The mantle’s isotopic compositions of ruthenium (delivered during Late Accretion) and molybdenum 
(delivered throughout accretion) lie along the meteoritic mixing line ^Dauphas et al 2004) l. This suggests that the 
isotopic nature of the projectiles did not change significantly before and after the Moon-forming impact. Given that 
^itoussi and Bourdon} ( |2012[ ) concluded, on the basis of silicon isotope compositions, that the enstatite chondrite 
contribution to the bulk Earth is < 15% using silicon isotope ratios, it is unlikely that the late accreted mass was 
dominated by enstatite chondrite projectiles. Moreover, Wang and Becker ( 2013| l also show a limited contribution of 
ordinary and enstatite chondrites to Late Accretion from an analysis of the Earth’s Se/lr and Te/Ir ratios. 

As stated at the beginning of this section, the number of bodies delivering most of Late Accretion must be small. 
Eor a single isotopic system, it is possible that the isotopic anomalies of the few largest projectiles relative to the 
Earth-Moon system cancel out. However, it is difficult for this to happen for two or more isotopic systems because the 
anomalies are different. Eor instance for ordinary chondrites and CR chondrites are opposite in sign and almost 

the same in absolute magnitude, but the eTi®° of CR chondrites is five times bigger than that of ordinary chondrites. 
Lunar-terrestrial tungsten anomaly: Given that the Earth and the Moon have identical (at analytic precision) isotope 

Spicuzza et al. 2007|l, Cr (Lugmair and 


compositions for many non-radiogenic elements (O {Wiechert et al. 2001 


Shukolyukov\\l9^^\Trinquier et al.^|2008|l, Ti ^Leya et al. 2008|l and Si (Armytage et al. 2012 Fitoussi and Bourdon 


20121), it is reasonable to assume that the Earth and the Moon started off with identical W compositions as well. 
Today, the difference in isotope composition between the Moon and the Earth is ( Touboul et al. 


20071 


,182, 


WMoon-e^“"WEarth= 0.09 ±0.1, 
where e^®^WMoon and e^®^WEarth are contributions to the tungsten anomaly after the Moon-forming impact until now. 


Each contribution can be broken into two categories: e^®^Wsampie = ± 


contributions 


from small and large bodies, respectively. We assume small, proportionately delivered bodies provide material of bulk 
chondritic composition. The small projectile contribution to the lunar tungsten anomaly is: 


('Hfl 

,182-i,rSmall ,182-i,r VW/Moon 
E Wmooh = e WchonyHfT- 
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= -4333 
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where e^®^Wchon = —2 is the average of the abundance of chondrites {Kleine et al. 20091, = 26 and 


= 1 are the halfnium/tungsten ratios of the Moon [Kleine et al. 20091 and the chondritic average iKleine 


^W/Chnn .. - r ■ - 

et a/.||2009i, respectively, me = 1.8 x 10 is the lunar late accreted mass [Day et al. 2007 Day and Walker\ 


20111 which is delivered entirely by small projectiles, and trim = 0.012 M® is the mass of the Moon’s mantle. 
Eor the Earth, the small projectile contribution to the terrestrial tungsten anomaly is: 


,182-ix;Small _ ,182 
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where (^)Earth “ Earth’s hafnium/tungsten ratio [Konig et al. 2011 1 , Ms = ame is the mass of small 

the 


projectiles impacting the Earth in terms of the lunar late accreted mass and the relative gravitational cross-section of 
the Earth and Moon since small projectiles are delivered proportionately, and Mm = 0.675 M 0 is the mass of the 
Earth’s mantle. 

Large bodies are stochastically delivered according to Bottke et al. (2010|l and thus presumably contributed only 


to the Earth. Thus they do not contribute to the Moon = 0. Earth contributions are complicated by 

the different possible fates of the projectile’s core, which is captured by the parameter X. If the projectile’s core is 
completely accreted in the Earth’s mantle X = 1, then the planetesimal contributes like a chondritic body; however if 
A < 1 then the tungsten anomaly contribution has two components: ± The hrst 

component corresponds to the fraction X of the projectile that merges with the Earth’s mantle delivering a chondritic 
contribution: 

tw)Earth f XMl \ 

V Mm J 
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where we used the relationship XMl = YMc — arric. 

The other component corresponds to the 1 — X fraction of each projectile’s core that sinks and merges with the 
Earth’s core. This component contributes only projectile mantle (achondritic) material to the terrestrial mantle tungsten 
anomaly; 


182\wAchon 
^ ^ Earth 


e'®'WAeho„ 

e^®"WAchon 

fHfj 

V W / Achon 


/ 0.675(1 -X)Ma \ 
f )ac.o„ I ) 

^ f {l-X) iYMc-amc) \ 

\ ) 


( 1 ) 

( 2 ) 


where we have assumed, for simplicity that the fraction of the projectile mantle to its total mass and the Earth’s mantle 
to its total mass are the same, i.e. 0.675. The formula above contains two unknown parameters for the mantle material 
of the large projectiles: and (w)Achon' analogs for the mantles of differentiated bodies come from 

the achondrite meteorite collections, and we list three examples in Table basaltic eucrites, aubrites and martian 
meteorites; 


Achondrites 





e'®^WAcho„ 

(Htj 

Basaltic Eucrites ( 

Kleine et al. 

2009 1 

22 

27 

Aubrites iPetitat et al. 

2008 

1 


11 

30 

Mars (Kleine et al 




0.4 

4.0 


Table 2: Tungsten abundances e^®^WAchon (column 2) and (^) Achon (column 3) for achondritic mantle materials 
(column 1 ). 


Eor each analog, we calculate the terrestrial tungsten anomaly contribution from differentiated mantles 
and then the tungsten difference between the Moon and the Earth e^®^WMoon — Wearth from the expressions above. 
This expression must match 0.09 ± 0.1 and is a function of X and Y. We overlay this tungsten difference constraint 
on top of the parameter space map shown in Panel a of Extended Data Eigure 4 to create new maps (see panels b-d) 
showing the contour that matches the nominal tungsten difference of 0.09 as a green solid line, as well as the l-cr limits 
as green dashed lines (regions exterior to the l-cr uncertainty are red). 

Regardless of the analog, the lunar-terrestrial tungsten difference places an upper and lower limit to the late ac¬ 
creted mass. It also always removes the possibility that the simple assumptions made in ^ottke ef q/.|(|2010|i ar e correct 
since the lunar-terrestrial tungsten difference is always too large when X = 1 and Y = 1. \Bottke et a(.| ( |2010 1 checked 
this constraint at the time of writing, but we are using here the updated Konig et al. (201 1| values for = 26. 

Using eucrite and the aubrite compositions rules out late accreted masses exceeding 0.013 and 0.033 Mq, respec¬ 
tively, but most of the parameter space consistent with the tungsten difference constraint implies a late accreted mass 
Mla ^ 0.01 M 0 , particularly when considering eucrite-like impactors. Only for a Mars-like composition is the con¬ 
straint looser. We stress that enstatite chondrites are very reduced, so a large enstatite body which underwent an early 
differentiation is expected to have a large hafnium-tungsten ratio and a large tungsten isotopic anomaly more similar 
to what is observed in the HEDs and aubrites. 


2.5 Conclusions on Late Accretion: 

The Earth mantle HSE budget, mantle tungsten isotope composition, the small difference in tungsten isotope compo¬ 
sition between the Moon and the Earth and other isotopic systems constrain the late accreted mass Ml a- Assumptions 
need to be made on the nature of the late accreted projectiles, but it turns out to be unlikely that the late accreted 
mass was larger than 0.01 M 0 . The standard model, which assumes that the late accreted mass is the same as the 
chondritic mass {Mla = Me = 4.8 ± 1.6 x 10“®), seems quite correct. Moreover it is also consistent with other 
non-geochemical constraints. In fact, Simone Marchi {personal communication) reports that in his Monte Carlo sim¬ 
ulations of the \Bottke et al.|p010| l scenario, calibrated on the lunar impact flux, the total mass delivered to the Earth 
results in less than 1% of the Earth’s mass in the vast majority of cases. Scenarios delivering more than 0.01 M 0 
in late accreted mass, while being potentially consistent with the constraints discussed in this paper, require ad-hoc 
assumptions on the chemical and isotopic properties of the dominant projectiles. 
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a) Mass-orbit distribution 


b) Time of last giant impact 



0.02 


0.5 1.0 2.0 

Semi-major axis (AU) 



c) Concentration statistic d) Angular momentum deficit 



Figure 3: Four panels comparing the classical (triangles) and Grand Tack (circles) scenarios. Panel a shows the mass of 
each planet in every simulated Solar System as a function of semi-major axis. Planets within grey zone are considered 
Earth-like. Planets at Mars’ distance from the Sun are too massive in classical simulations but are correct in Grand 
Tack simulations. The next three panels show the relative Late Accretion mass as a function of different parameters 
and the chondritic mass and uncertainty Me =4.8 ± 1.6 x 10“^ M 0 as a dashed line and grey region. Panel b is the 
same as Figure 1 of the main text, but with a broader time scale. Panel c and d show the relative Late Accretion mass 
and two orbital structure statistics; concentration Sr. and angular momentum deficit Sd- The concentration statistic 

(: ^ 


is defined [Chambers 


20011; Sc = max | 


^^ T,j'mj{Logio{a/aj))^ 

and Gj are the mass and semi-major axis of the jth planet. Sc 


, where the summation is over the planets and rrij 
is then the maximum of this function as a is varied. 


The Solar System has Sc = 89.9 and in panel c is marked with a vertical dashed line. Grand Tack simulations 
are as or more concent rated than classical simulations. The angular momentum deficit statistic is defined ^Laskar\ 

Wil: Sd = 


1997 


Chambers 


J2j rrij^/a- 


-, where the summation is again over the planets and Cj and 


are the eccenfncity and inclination of the jth planet. The Solar System has Sd — 0.0018 and this marked by a 
vertical dashed line in panel d. A late giant planet instability (^500 My after first Solids in the Solar System) likely 
excites the inner Solar System, since these simulations end before the instability they should only be a fraction of the 
current value ^Brasser et a/.|[20T3| l. Grand Tack simulations are either as dynamically excited or colder than classical 
simulations, although they are in some cases hotter than the current terrestrial planets. 
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a) Total embiyo to planetesimal mass ratio 1:1 b) Total embryo to planetesimal mass ratio 2:1 



Time of last giant impact (Myr) Time of last giant impact (Myr) 


c) Total embiyo to planetesimal mass ratio 4:1 d) Total embiyo to planetesimal mass ratio 8:1 



Time of last giant impact (Myr) Time of last giant impact (Myr) 


Figure 4: Figure 1 is reproduced in each panel, showing the Earth-like planet relative late accreted mass as a function 
of last giant impact time for each scenario within the classical category. For reference, the horizontal dashed line and 
enclosing darkened region are the best estimate and 1 -ct uncertainty of the late veneer mass Mq = 4.8 ± 1.6x 10“^ 
M 0 . Each panel shows the Earth-like planet relative late accreted mass as a function of last giant impact time for 
different initial ratios of the total mass in embryos to planetesimals including 1/1 in a, 2/1 in b, 4/1 in c, and 8/1 in 
d. The mean time of the last giant impact increases and the relative Late Accretion mass decreases as the initial total 
embryo to planetesimal mass ratio increases. 
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Mantle accretion fraction of projectile cores, X 
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Figure 5; The late accreted mass Mla shown as grayscale contours of the two parameters; X and Y. The blue dashed 
contours indicate the location of a few specific late accreted masses and are labeled in the framed boxes. For this 
figure, we assumed that the Earth receives 52 times as much mass in small projectiles than the Moon (cr = 52). The 
red region is inaccessible, since the measured chondritic mass in the lunar mantle requires a minimum flux of small 
planetesimals onto the Earth’s mantle. 


18 















a) Accretion history 


b) Basaltic eucrites 
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Figure 6; Each panel is like Extended Data Eigure 3 but also showing additional constraints. The constraint X > 0.53, 
deduced by a reanalysis of Rudge et al.\^20l^ with the Hf/W ratio updated from Konig et al. (20111, is shown in 
light red in each panel. In panels b-d, the solid green contours indicate the set of X, Y parameters that reproduce 


the nominal difference in W-isotope composition between the Moon and the Earth (0.09 in e units; see Touboul et al 


( |2007| l). The difference among these panels is in the assumed Hf/W and W-isotope composition for the mantle of the 
differentiated projectiles; panel b assumes values typical of eucrites; panel c is for aubrites and panel d is for Mars. 
The green dashed contours are the 1-cr uncertainties on X, Y related to the uncertainty on the difference between the 
W-isotope composition between the Moon and the earth: zbO.l in e units. The areas exterior to these 1 — a uncertainties 
have been colored red. 
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